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Abstract 

We provide here an alternative derivation of the generalization of the nonlinear Turin model for dispersion 
"^q". unmanaged coherent optical links provided in Johannisson's report (T). 
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r/3 1 I. INTRODUCTION 

. 1 Goal of this paper is to provide a simplified derivation of the results appearing in the recent ArXiv posting of P. 
Johannisson (H on a generalization of the well-known frequency-domain nonlinear interference (NLI) analytical 
! model for dispersion unmanaged (DU) coherent systems introduced by Turin's group in O. 
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II. Frequency Domain NLI RP1 Solution 

1 We start from the dual-polarization (DP) single-channel first-order Regular Perturbation (RP1) solution of the 
ON ■ dispersion-managed nonlinear Schroedinger equation (DMNLSE) ( [3], Appendix 2): 

o 



X 



where L is the total link length, and the NLI perturbation field is 

u p (l, f) = -jp ir ic(hf 2 )u(o, f + fi)ir\o, f + h + h)u(o, f + f 2 )f 1 h (2) 



where: 

i) boldface fields are 2x1 vectors containing the Fourier transforms (denoted by a tilde) of the X and Y polar- 
izations in the transmitter polarization frame of reference; a dagger stands for transposition and conjugation; and 
the DP field power is normalized to an arbitrary reference power Pq (which in (T| is chosen as the per-polarization 
average power); 

ii) the un-normalized scalar frequency kernel is defined as (Cfr. 12 eq. (45) and Q eq. (25)): 

K(F)± [%'( s )g(s)e-^ s ^ 2F ds (3) 
Jo 

where F = f\f 2 is the product of two frequencies, 7' = I7 with 7 the fiber nonlinear coefficient, Q(s) the line 
power gain from z = to z = s, and C(s) = — f3 2 (z)dz is the cumulated dispersion in the transmission fibers 
(with dispersion coefficient /3 2 ) up to coordinate s. We choose here the frequency-normalizing rate in [|3] as R = 1, 
i.e., we do not normalize the frequency axis, as in HI. In Q we use the normalized kernel 

l{ ) /C(0) 
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which then at F = equals 1 . The nonlinear phase referred to power Pq is 

<5> NL = P K(p). 
Hence by multiplying and dividing by /C(0) we can recast (O as 

U P (L, f) = -j$ NL If™ f,{hf 2 )U{Q, f + h)u\0, f + fi + f 2 )U(0, f + f 2 )iih (4) 
From the X component of the RP1 solution writes explicitly as 

Ux ' p( ^ J) = tr f,(hf 2 )u x (o, f + A)u* x (o, f + h + h)u x (o, f + / 2 )f 1 f 2 + 

-J®NL JJ-oo 

oo 

ti(flf2)Ux(P, f + /l)^(0, f + fl + f2)Uy(0, f + f2)ht 2 (5) 



-oo 

where the first line gives the self -phase modulation (SPM) of X on X, while the second line gives the intra-channel 
cross-polarization modulation (I-XPolM) of Y on X. A perfectly dual expression for component Y is obtained by 
exchanging the indices x and y. 

III. Gaussian Assumption and Johannisson's Result 

In CD, HI the key assumption is that the input fields are composed of independent spectral lines with Gaussian 
amplitude: 



U x (0, f) = VTo Yl MG x (kfo)8(f - kf 



k=—oo 



U y (0, /) = <//££ (k\JG y (kf )8(f - kfo) 

A:= — OO 

with £k and independent identically distributed standard (i.e. zero-mean unit variance) circular complex Gaus- 
sian random variables (RV). Such signals do have a per-polarization power spectral density G x / y {f) (normalized 
to Po) in the limit fo — > Q. Then after long statistical averaging calculations, one gets the power spectral density 
of the U X:P (L, f) RV as ( d, eq. (89). Note that our PSD G{f) is normalized such that G(f) = P G(f), where G 
is the un-normalized PSD per polarization. Also, P x = Po G x (f)df and P y = Po G y {f)df) 



G x , p (f) = p o{2 // \mh - f){h ~ fW G x (h)G x (f 2 )G x (h + f 2 - f)^ 
|iC((/i - f)(f 2 - f))\ 2 G x {h)G y {f 2 )G y {h + f 2 - f)U 2 

-oo 

(roo roo roo roo \ 

4( / G x {f)dff + 4 / G x {f)df / G y {f)df + ( / G y (f)df) 2 } (6) 
J —oo J — oo J —oo J —oo J 

and a dual expression for Y is obtained by swapping x -B- y. Recall that G XjP (f) is the NLI PSD, normalized by 
Po- 

An equivalent form of © is the following 

"OO 

|2 



G x , p (f) =Hl{2 // WihW G x (f + h)G x {f + f 2 )G x (f + f 1 + f 2 )U 2 

JJ — OO 

rroo 

+ // \v(hf2)\ 2 G x (f + h)G y {f + / 2 )G y (/ + / 1 + / 2 )f 1 f 2 

(roo roo roo roo \ 

4(/ G x (f)dff + A G x {f)df Gy (f)df + ( G y (f)df) 2 )} (7) 
J —oo J —oo J —oo J —oo / 

which better shows the formal parallel with the field equation (f5]): the field double integral in fi, f 2 of the product 
kernel-field-field* -field becomes a PSD double integral in /i, f 2 of the product squared kernel magnitude-PSD- 
PSD-PSD. 

It is the purpose of the remaining part of this paper to provide a new proof of (0. 
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IV. The New Proof 

We now start from (f5]) and make the following two assumptions regarding the input X,Y fields U x (0, t), U y (0, t): 

1) they are wide-sense stationary (WSS); 

2) they are jointly Gaussian processes. 

Regarding assumption 1), we plan to exploit the following extension of result ( (U, p. 418, eq. (12-76) ): 
Theorem 1 

Consider the jointly WSS stochastic processes x(t) and y(t), and let 

/oo 
-oo 

the Fourier transform of x ( in the mean-square (MS) sense), and Y(f) is similarly defined. Let their cross power 
spectral density (PSD) be G xy (f) = F[R xy {T)] = F[E[x{t + r)y*(t)]]. Then 

E[X{f)Y*(u)} = G xy (f)5(f -u) = G xy {f)5{u - f) □ (8) 

As a byproduct, we also have 

E[X(f)X*(v)] = G x (f)S(f -u) = G x (f)5(u - f). 

This theorem thus shows that the Fourier transform of any MS-integrable WSS process is nonstationary white 
noise, and thus the spectral lines of its Fourier transform are uncorrelated. 

Regarding assumption 2), we plan to exploit the following result, known as the complex Gaussian moment 
theorem (CGMT), a generalization to complex variables of Isserlis theorem 151, j6ll: 

Theorem 2 

Let U\, U2, ■ ■ ■ , be zero-mean jointly circular complex Gaussian random variables. Then 

E[UlUl..UtU k+1 U k+2 ...U 2k ] = £ E[UlU p ]E[U$U q \ E[U* k U r ] (9) 

where denotes a summation over the k\ possible permutations (p, q, r) of indices (k+1, k+2, 2k) □ 
For instance, 

E[UlU^UiU 4 U 5 U 6 ] = E[UZU 4 ]E[UtU 5 ]E[U£U 6 ] 
+ E[U^U 4 ]E[U^U 6 ]E[U^U 5 ] 
+ E[U^U 5 ]E[U^U A ]E[UIU 6 ] 
+ E[UtU 5 ]E[UZ U 6 ]E[UIU 4 ] 
+ E[U^U 6 ]E[U^U 4 ]E[U^U 5 ] 

+ E[UIU 6 ]E[U^U 5 ]E[U^U 4 ]. (10) 



Let's now start the new proof. We are interested in the PSD G XtP (f) of the NLI field U XtP (L, t) = [U X)P (L, /)]. 
By theorem 1 we have: 

E[U X , P (L, f)U* p (L,u)] = G x>p (f)S(u - f). (11) 
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The left hand side can be explicitly calculated using (T5]): 

E[U x>p (L,f)U* tP (L,u)] 

rrco 

E[ // ^(/i/ 2 )[C4(0, / + /i)E£(0, / + /i + / 2 )C4(0, / + / 2 )+ 

^(0, / + fl)U;(0, f + fl + h)Uy(0, f + / 2 )]f!f 3 - 

vihhTP* (0, u + f 3 )U x {0, u + h + f 4 )U* (0, u + / 4 )+ 
C/*(0, « + h)U y {0, u + h + h)U;(0, u + / 4 )]f 3 f 4 ] = 

i 

{E[U x (0, f + /i)e£(0, / + h + f 2 )U x (0, f + / 2 )t£(0, u + / 3 )C4(0, u + / 3 + / 4 )^(o, « + / 4 )]+ 

£[e>*(o, / + h)u;(o, f + h + f2)u v (p, f + h)u*(o, u + h)u y (o, u + f 3 + h)u;(o, u + h)]+ 

2Re(£[^(0, / + fi)U*(0, f + fi + f 2 )U x (0, f + f 2 )U* x (0, u + f 3 )U y (0, u + f 3 + / 4 )E£(0, u + / 4 )])}. (12) 

Now, putting together Theorems 1 and 2, Appendix 1 shows the following 
Theorem 3 

For jointly stationary circular complex Gaussian zero-mean processes A(t), B(t),C(t), D(t), E(t), F(t) we 
have the general formula 

' Mf + fi)B*(f + fi + f 2 )C(f + f 2 )D*(u + h)E(u + h + U)F*{u + U) 

G ab (f + h)G cd {f)G ef {f + h)5{h)8{h)+ 

Gabif + h)G ed {f + h)G cf (f)5(f 2 )S(U)+ 
G cb (f + h)G ad {f)G ef {f + U)5{h)5{h)+ 
G cb (f + f 2 )G ed (f + h)G af (f)5{h)5{U)+ 
G eb (f + + h)G ad (f + h)G cf (f + f 2 )5(f 3 - fi)5(h ~ h)+ 

Gebif + fl + h)GalU + f 2 )G af {f + fl)5(U ~ fl)5(f 3 - f 2 ) 

■8{u -f) □ 

We next apply the general formula ( fT3l to the three expectations in (fT2l to get: 
First expectation: 

E [U x (0,f + fi)U*(0J + fi + f 2 )U x (0J + f 2 )U* x (0,u + f 3 )U x (0, u + f 3 + f 4 )U*(0,u + h) 

s(u - f) [c xx (f + j\ )(;.,., if s;.,., if + h)5(f 2 )s(f 3 )+ 

G xx (f + h)G xx (f + h)G xx (f)5(f 2 )5(U)+ 
G xx (f + f >)(:.,., if )(!,■, -if + h)S(fi)5(f 3 )+ 
G xx (f + f 2 )G xx (f + f 3 )G xx (f)8(f 1 )8(f 4 )+ 
G xx (f + + h)G xx (f + h)G xx (f + f 2 )5(f 3 - h)5{f 4 ~ h)+ 



(13) 



G xx (f + h + f 2 )G xx (f + f 2 )G xx {f + h)5{U - h)8{f 3 - f 2 ) 



(14) 
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where G xx = G x . Second expectation: 



E 



U x (0, f + /i)C/*(0, f + h + f 2 )U y (0, f + f 2 )U*(0, u + f 3 )U y (0, u + h + h)U;(0, u + / 4 ) 
S(u - f) [G xy U + h)G yx (f)G yy (f + f 4 )5(h)S(h)+ 
G xy (f + h)G yx {f + h)G yy (f)5(f 2 )6(U)+ 

G yy (f + f 2 )G xx (f)Gyy(f + U)5(fl)S(f 3 ) + 

G yy (f + f 2 )G yx (f + h)G xy (f)6(f 1 )S(f 4 )+ 
G yy (f + /i + / 2 )G,,(/ + /i)G w (/ + / 2 )5(/ 3 - h)5(U - h)+ 
Gyy(f + h + h)G yx {f + f 2 )G xy (f + h)5{h - fi)S(f 3 - h) 

where G yy = G y , and assuming uncorrected X and Y we get 



(15) 



E 



Ux(0, f + /i)CC(0, f + fi + h)U y (0, f + f 2 )U*(0, u + f 3 )U y (0, u + f 3 + / 4 )tL*(0, n + / 4 ) 



5(« - /) [G w (/ + f 2 )G xx (f)G yy (f + U)S(h)6(f 3 )+ 
G yy (f + h + /2)G M (/ + /i)G w (/ + / 2 )«5(/ 3 - /i)<5(/ 4 - / 2 ) 
Third expectation: 

£ [o*(o, / + fi)u*(o, f + h + h)u x (o, f + h)u* x (o, u + f 3 )u y (o, u + f 3 + U)u;(o, u + u) 

6(u - f) [G xx (f + h)G xx {f)G yy (f + U)5{f 2 )8(h)+ 

G xx {f + h)G yx (f + h)G xy (f)8(f 2 )5(U)+ 
G xx (f + f 2 )G xx (f)G yy (f + f A )5(f l )5(h)+ 
G xx (f + f 2 )G yx (f + f 3 )G xy (f)5(f 1 )8(h)+ 
G yx (f + h + h)G xx (f + h)G xy {f + f 2 )5(f 3 - h)5(U - h)+ 
G yx {f + + h)G xx {f + f 2 )G xy {f + h)5(U - h)5{h ~ ft) 
and assuming uncorrected X and Y we get 

E [U x (0, f + fi)U*(0, f + fi + h)U x (0, f + h)U* x (0, u + h)U y (0, u + f 3 + / 4 )tf*(0, u + / 4 ) 
S(u - f) [G xx (f + h)G xx {f)G yy (f + h)S(f 2 )S(h)+ 
G xx (f + f 2 )G xx (f)G yy (f + U)S(h)S(f 3 ) 
Substitution of (|T4]l,(|T6l),(|T8l» into (O finally gives 
E[U x , p (L,f)U* jP (L,u 



(16) 



(17) 



(18) 



$2 



S(u-f){ 



flf 2 f3f^(/l/2)r?(/3/4)*- 



■ { [G x (f + h)G x (f)G x (f + h)S(f 2 )S(f 3 ) + G x (f + h)G x {f + f 3 )G x (f)5(f 2 )5(h)+ 
G x (f + f 2 )G x (f)G x (f + h)5(ft)S(f 3 ) + G x (f + f 2 )G x (f + h)G x {f)5{h)5(h)+ 
G x {f + fi + f 2 )G x (f + h)G x {f + f 2 )5{h - fx)S{U - f 2 )+ 



+ 



G x (f + + f 2 )G x (f + f 2 )G x (f + h)5{U - hWz ~ h) 

~G y (f + f 2 )G x (f)G y (f + h)S(fiWs) + G y (f + + f 2 )G x (f + h)G y {f + f 2 )5(f 3 - fi)5(h ~ h) 
2Re( \G x (f + fi)G x (f)G y (f + U)5(f 2 )5(f 3 ) + G x (f + f 2 )G x (f)G y (f + U)S(fi)S(f 3 )] )}}. 



+ 
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From (TTTT >. the term multiplying S(u — f) must be the desired PSD. Each pair of delta removes two integrals, so 
that the PSD turns out to be (first two lines above produce first 4 lines, third line above produces 5th line, 4th line 
above produces 6th and 7th lines, and final line above produces the last two lines): 



oo 

} =l 7 ?(°)l 2 / G x (f + h)G x (f)G x (f + f^f, 

— oo 
oo 

+ l*K°)l 2 / G x (f + h)G x (f + f 3 )G x (f)tits 



$2 

NL 

— oo 

oo 



— oo 
oo 



+ \fi(0)\ 2 [f G x (f + f 2 )G x {f)G x {f + h)t 2 U 

— oo 
oo 

+ l^(°)l 2 / G M + + h)GM)hh 

— oo 

oo 

+ 2 ff mif 2 )\ 2 G x (f + f 1 + f 2 )G x (f + h)G x {f + f 2 )f 1 ? 2 



-oo 

'DO 



+ I^O)! 2 ^ Gy(f + f 2 )G X (f)Gy(f + / 4 )f 2 f 4 

— oo 

oo 

+ jj mif 2 )\ 2 Gy(f + A+ / 2 )4(/ + h)Gy{f + h)i A i 2 

— OO 

(OO 
l^(0)| 2 jj G x (f + h)G x {f)G y U + / 4 )f 1 f 4 + 
— oo 

oo \ 

|^(0)| 2 J G , ,(/ + / 2 )G.(/)G ?/ (/ + /4)f 2 f4 • 

— oo / 

In summary, considering that by construction 77(0) = 1, we have: 

00 

%^ = 2 ft Wih)\ 2 G x {f + f l + f 2 )G x U + h)G x {f + f 2 )fj 2 

—00 
00 

+ jj mif2)\ 2 Gy(f + f 1+ f 2 )G X (f + h)Gy{f + f 2 )f A f 2 

— OO 

(f'QO POO POO POO \ 

4( / Gx(/)#) 2 + 4 / G x (f)df / Gj,(/)d/ + 1( / G,(/)(i/) 2 
J — 00 i/ — 00 i/ — 00 i/ — 00 / 

which confirms Johannisson's equation © and completes the desired alternative proof. 



V. Conclusions 

We have presented an alternative derivation of Johannisson's result JT|. We first remark that our new method is 
able to deal with correlated X and Y, although this feature was not exploited in the present paper. Next we note that 
we did not have to assume independent input spectral lines: this comes naturally from the stationarity of the input 
process. Finally, the truly critical assumption in the model in (H, Q is therefore the assumption of Gaussianity 
at any z during propagation, which is implicit in the assumption of a Gaussian input process, and the fact that the 
"forcing terms" in the RP equation are the linearly distorted signals at any z, which thus remain Gaussian. 



7 



Therefore the true limit of the model in 12, (21 is that indeed starting from a non-Gaussian spectrum such 
as the one of a digitally modulated signafl it takes some finite propagation in a non-infinite dispersion line to 
approximately get both a Gaussian spectrum and a Gaussian-like time-domain signal. 

Appendix 1 

In this Appendix we prove Theorem 3 in the text. Assuming jointly stationary circular complex Gaussian zero- 
mean processes A(t),B(t), C(t),D(t), E(t),F(t), we have by using © in ( fTUb : 

T^E [A(f + fx)B*{f + h + f 2 )C(f + f 2 )D*(u + h)E(u + f 3 + U)F*(u + / 4 ) 

E[B*(f + h + f 2 )A(f + h)] E[D*(u + h)C(f + f 2 )\ E[F*(u + U)E{u + / 3 + / 4 )]+ 

'^•^^^^^^^^^^^^^^K. s-^^^^^^^^^^^^^^—^ S-^^^^^^^^^^^S ^^^^^^^^^^^^^^^^ J 

G ab (f+h)S(f 2 ) G ad {f+f 2 )S(u+f 3 -f-f 2 ) G e/ (u+/ 3 +/ 4 )<5(/ 3 ) 

E[B*(f + h + f 2 )A(f + /i)] E[D*(u + h)E{u + f 3 + / 4 )] E[F*(u + f 4 )C(f + f 2 )} + 
> , ' v v ' v , ' 

G ab (f+fi)S(f 2 ) G ed (u+f s +U)S(U) Gc/ (/+/ 2 )5(«+/ 4 -/-/ 2 ) 

E[B*(f + h + / 2 )C(/ + / 2 )] + / 3 )i(/ + /i)] + / 4 )^(« + /a + / 4 )] + 

V v ' V v ' V v ' 

G ab (f+f 2 )S(h) G ad U+h)5{u+h~f-fi) G ef (u+f 3 +U)S(f 3 ) 

E[B*(f + h + f 2 )C(f + f 2 )} E[D*(u + h)E{u + h + U)} E[F*(u + f 4 )A(f + 

G ci) (/+/ 2 )5(/ 1 ) G Ed («+/ 3 +/ 4 )<5(/ 4 ) G a/ (/+/i)<5(«+/4-/~/i) 

E[B*(f + fi + f 2 )E(u + h + h)] E[D*(u + f 3 )A(f + ft] E[F*(u + f 4 )C(f + f 2 )] + 

V v ' V v ' V v ' 

G e6 (u+/ 3 +/4)a(/+/ 1 +/ 2 -u-/ 3 -/4) G ad (/+/!)5(«+/ 3 -/-/i) G cf (f+f 2 )S(u+U-f-f 2 ) 
E[B*(f + h + f 2 )E(u + f 3 + / 4 )] + / 3 )C(/ + / 2 )] £[F*(u + / 4 )i(/ + /,)] = 



G e6 («+/ 3 +/ 4 )5(/+/i+/2-«-/3-/ 4 ) G c<j (/+/ 2 )<5(m+/ 3 -/-/ 2 ) G o/ (/+/i)5(«+/4-/-/0 



thus 



r = G a6 (/ + /ijGrfC/ + / 2 )G e/ (u + / 3 + h)6(h)6(f 3 )6(u + f 3 -f-f 2 ) 
+ G ab (f + /i)Ged(« + /s + /4)G c/ (/ + f 2 )8(f 2 )S(U)S(u + U-f-h) 
+ G cb (f + f 2 )G ad (f + h)G ef (u + h + U)S(fi)S(f 3 )S( U + /3-/-/1) 
+ <?<*(/ + / 2 )G ed (u + / 3 + / 4 )G a/ (/ + fi)6(fi)5(U)6(u + h-f-h) 
+ G eb (u + h + U)G ad (f + fi)G cf (f + f 2 )- 
■S(u + h-f- f 2 )5(u + h-f- h)8(J + h + f 2 -u-h-U) 
+ G eb (u + h + h)G cd (f + h)G a fU + h)- 

■ 8{u + h-f- h)8{u + h-f- f 2 )S(f + f l + f 2 - u -h- U). 
Now we use the sampling property of the delta to write, e.g. for the first line where f 2 = and / 3 = 0, 

Gab(f + fi)G cd (f)G ef (u + h)5(f 2 )5(h)5{u - f) 
and e.g. for the last line where u + / 4 = / + f\ and u + / 3 = / + f 2 which we add up to get 

u + / 3 + U = (f ~ u) + f + fi + f 2 

whence 

f+h+h-u-h-h=u-f 

so that the last line writes as 

G eb {u + h + f 4 )G cd (f + h)G af (f + fi)S(u + U-f- h)S(u + h-f- f 2 )5(f + f 1 + f 2 - u -h-U) = 

G eb ((f -u) + f + h + h)G cd {f + f 2 )G af (f + h) .5(u + h-f- h)5(u + f 3 -f- f 2 )5(u - f) (US£ = =/) 
G eb (f + fi + h)G cd (f + f 2 )G af (f + h)5{U - fi)S(f s - f 2 )5(u - f). 

'Although the authors in J2) present in their Appendix B an appealing heuristic justification of their Gaussian signal assumption, still 
their invoking the central limit theorem at their equation (37) is not rigorous. They would conclude that any digitally modulated signal with 
any number of levels has a Gaussian Fourier transform (which in turn implies the time-domain signal itself is Gaussian), which is clearly 
not the case. 



We therefore get 

T = G ab {f + h)G cd (f)G ef (f + h)5{f 2 )5{h)5{u - f) 
+ G ab (f + h)G ed (f + h)G c f(f)5{f 2 )5{h)5{u - f) 
+ G cb (f + f 2 )G ad (f)G ef (f + h)S(f 1 )8(fs)S(u - f) 
+ G cb (f + h)Ged{f + h)G af {f)8{h)5{U)8{u - f) 
+ G eb (f + f l + h)G ad {f + h)G cf U + /2W/4 - f2)S(h - h)S(u - f) 
+ G eb (f + + f 2 )G cd (f + h)G a} {f + /i)«5(/ 4 - /i)<5(/ 3 - f2)5(u - f). 

whence the final form (fT3l) given in Theorem 3. 



